Accepted by ApJ 

Preprint typeset using I^T^X style cmulatcapj v. 08/22/09 



THE SPITZER C2D SURVEY OF NEARBY DENSE CORES. IX. DISCOVERY OF A VERY LOW 
LUMINOSITY OBJECT DRIVING A MOLECULAR OUTFLOW IN THE DENSE CORE L673-7 

Michael M. Dunham 1,2 , Neal J. Evans II 1 , Tyler L. Bourke 3 , Philip C. Myers 3 , Tracy L. Huard 4 , & Amelia 

M. Stutz 5,6 

Accepted by ApJ 

ABSTRACT 

We present new infrared, submillimeter, and millimeter observations of the dense core L673-7 and 
report the discovery of a low-luminosity, embedded Class protostar driving a molecular outflow. 
L673-7 is seen in absorption against the mid-infrared background in 5.8, 8, and 24 /zm Spitzer images, 
allowing for a derivation of the column density profile and total enclosed mass of L673-7, independent 
of dust temperature assumptions. Estimates of the core mass from these absorption profiles range 
from 0.2 — 4.5 M©. Millimeter continuum emission indicates a mass of ~ 2 M©, both from a direct 
calculation assuming isothermal dust and from dust radiative transfer models constrained by the 
millimeter observations. We use dust radiative transfer models to constrain the internal luminosity of 
L673-7, defined to be the luminosity of the central source and excluding the luminosity from external 
heating, to be Li nt = 0.01 — 0.045 L©, with Li n t ~ 0.04 L© the most likely value. L673-7 is thus 
classified as a very low luminosity object (VeLLO), and is among the lowest luminosity VeLLOs yet 
studied. We calculate the kinematic and dynamic properties of the molecular outflow in the standard 
manner. From the outflow properties and standard assumptions regarding the driving of outflows, 
we calculate the time-averaged protostellar mass accretion rate, total protostellar mass accreted, and 
expected accretion luminosity to be (M acc ) > 1.2 x 10~ 6 ^^ M© yr _1 , Af acc > 0.07^ M© and 
Lace > 0.36 L©, respectively. The discrepancy between this calculated L acc and the Li n t derived from 
dust radiative transfer models indicates that the current accretion rate is much lower than the average 
rate over the lifetime of the outflow. Although the protostar embedded within L673-7 is consistent 
with currently being substellar, it is unlikely to remain as such given the substantial mass reservoir 
remaining in the core. 

Subject headings: ISM: individual (L673-7) - stars: formation - stars: low-mass, brown dwarfs 



1. INTRODUCTION 

The Spitzer Space Telescope Legacy Project "From 
Molecular Cores to Planet Forming Disks" (c2d; Evans 
ct al. 2003) has completed a 3.6 — 160 fim imaging sur- 
vey of nearby, low-mass star-forming regions, including 
approximately 100 small, dense starless and protostel- 
lar cores. The very first starless core observed by c2d, 
L1014, was found to harbor an embedded protostar with 
Lint ~ 0.09 L© (Young et al. 2004), where L int , the inter- 
nal luminosity, is the luminosity arising from the central 
protostar and disk and excludes luminosity from exter- 
nal heating by the interstellar radiation field (ISRF). The 
discovery of L1014 led to the definition of very low lu- 
minosity objects (VeLLOs; Di Francesco et al. 2007) as 
objects embedded within dense cores with Li nt < 0.1 L©. 

L1014 is not the only VeLLO discovered by c2d. Dun- 
ham et al. (2008) identified 15 candidate VeLLOs in the 
full c2d sample. In addition to L1014, three others have 
been studied in detail: IRAM 04191+1522 (hereafter 
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IRAM04191) 7 {L int ~ 0.08 L©; Dunham et al. 2006), 
L1521F {L mt ~ 0.05 L©; Bourke et al. 2006; Terebey et 
al. 2009), and L328 (L int ~ 0.05 L©; Lee et al. 2009). 
Other groups have also found objects consistent with 
the definition of VeLLOs: Chamaeleon-MMSl, either a 
VeLLO with Li nt ~ 0.01 — 0.02 L© or possibly a first hy- 
drostatic core (FHSC) (Belloche et al. 2006), SSTB213 
J041757, a candidate embedded proto brown dwarf with 
L int ~ 0.003 L© (Barrado et al. 2009), and L1448 IRS2E, 
a candidate FHSC with L bot < 0.1 L© (Chen et al. 2010). 
Assuming spherical mass accretion at the rate predicted 
by the standard model (M acc - 2 x 10~ 6 M© yr" 1 ; Shu, 
Adams, & Lizano 1987) onto an object with a typical 
protostellar radius of R ~ 3 R©, a protostar located on 
the stellar/substellar boundary (M — 0.08 M©) would 
have an accretion luminosity, L acc = GMM acc / R, of 
L 1.6 L©. VeLLOs, with luminosities more than an 
order of magnitude lower, must either feature mass accre- 
tion rates lower than predicted by the standard model, 
masses below the stellar/substellar boundary, or some 
combination of the two. 

Despite their similar luminosities, the VeLLOs stud- 
ied to date show very different core and molecular out- 
flow properties. L1014, IRAM04191, and L1521F all 
have substantial mass reservoirs remaining in their en- 

7 Unlike most VeLLOs, IRAM 04191+1522 was previously 
known to exist based on the detection of warm dust emission and 
a molecular outflow (Andre et al. 1999). 
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velopes (M env ^2 — 5 M Q ; Young et al. 2004; Dunham 
et al. 2006; Bourke et al. 2006), whereas L328 appears 
to only have M env ~ 0.1 M Q . IRAM04191 drives an 
extended, bright, bipolar molecular outflow detected in 
single-dish observations (Andre et al. 1999). L1014 drives 
a weak, compact outflow detected by the Submillimeter 
Array (Bourke et al. 2005) but not by single-dish ob- 
servations (Crapsi et al. 2005). L1521F and L328 both 
show some hints of outflows in their infrared morpholo- 
gies (Bourke et al. 2006; Lee et al. 2009), but neither 
show unambiguous evidence for driving outflows. Addi- 
tional VeLLOs must be studied in detail to learn more 
about the properties of this class of objects. 

VeLLOs represent one extreme of the protostellar lu- 
minosity distribution, which recent studies have shown is 
strongly peaked at low luminosities (Dunham et al. 2008; 
Enoch et al. 2009a; Evans et al. 2009). These results 
represent a confirmation and re-statement of the clas- 
sic "luminosity problem" whereby accretion at the stan- 
dard rate produces accretion luminosities higher than 
typically observed for embedded protostars (Kenyon et 
al. 1990). The luminosity distribution is inconsistent 
with a constant protostellar mass accretion rate through- 
out the duration of the protostellar phase and instead 
suggests that it is highly time variable with extended 
periods of very low accretion, as first suggested by 
Kenyon & Hartmann (1995). Vorobyov & Basu (2005a, 
2006, 2009) presented simulations showing that mate- 
rial will pile up in a circumstellar disk until the disk be- 
comes gravitationally unstable and dumps its mass onto 
the protostar in short-lived accretion bursts. Vorobyov 
(2009) compared the distribution of mass accretion rates 
in these simulations to those inferred from the luminosi- 
ties of protostars in three molecular clouds compiled by 
Enoch et al. (2009a) and concluded that their simulations 
reproduced some of the basic features of the observed 
distribution of mass accretion rates. Based on these sim- 
ulations, Dunham et al. (2010) presented evolutionary 
models modified to include a simple treatment of episodic 
mass accretion in the standard model of the collapse of a 
singular isothermal sphere (Shu 1977; Shu, Adams, & 
Lizano 1987; Young & Evans 2005) and showed that 
models with episodic accretion mostly resolved the lu- 
minosity problem and improved the match between the 
observed and predicted protostellar luminosity distribu- 
tions. 

Some direct observational evidence for episodic mass 
accretion in the embedded phase exists in the form of 
accretion bursts in Class I sources (e.g., Acosta-Pulido 
et al. 2007; Kospal et al. 2007; Fedele et al. 2007), a mis- 
match between the measured accretion rates onto the 
disk and protostar of NGC 1333-IRAS 4B (Watson et 
al. 2007), and episodicity in jets ejected from some pro- 
tostellar systems (e.g., Lee et al. 2007). Additionally, 
evidence for non-steady mass accretion can be found by 
studying the molecular outflows driven by protostars. 
Since outflows are driven by accretion they retain in- 
formation on the accretion history. The outflow prop- 
erties can be used to calculate time-averaged mass ac- 
cretion rates and thus accretion luminosities that can 
be directly compared to the protostellar internal lumi- 
nosities. Two Class sources studied in this manner, 
IRAM04191 and L1251A-IRS3, both have accretion lu- 
minosities inferred from their outflows much greater than 



their internal luminosities, possibly indicating lower cur- 
rent accretion rates compared to the time-average (Andre 
et al. 1999; Dunham et al. 2006; Lee et al. 2010). 

In this paper we present new infrared, submillimeter, 
and millimeter observations of the dense core L673-7 and 
report the discovery of a low-luminosity, embedded pro- 
tostar driving a molecular outflow. We provide a brief 
introduction to L673-7 in fj2j A summary of the obser- 
vations and data reduction is given in $3l We present a 
qualitative discussion of the results in SJU including the 
discovery of the protostar and the detection of L673-7 in 
absorption against the mid-infrared background in §4.11 
a discussion of the evolutionary status of the protostar 
in H4.2[ a discussion of the L673-7 core mass in 34.31 and 
the discovery of the molecular outflow in 34.41 In ij5]we 
use the absorption profile in the mid-infrared to derive 
the column density profile and total enclosed mass of 
L673-7. Dust radiative transfer models used to constrain 
the protostellar internal luminosity and confirm its sta- 
tus as a VeLLO are discussed in 01 In §T[ we calculate 
the kinematic and dynamic properties of the outflow and 
use these properties to calculate the timc-avcraged mass 
accretion rate onto the protostar and expected accretion 
luminosity from L673-7. Finally, SJU discusses the results 
and SJH presents our conclusions. 

2. L673-7 

L673-7 is part of the L673 cloud complex from the 
Lynds (1962) catalog, which lists L673 as an Opacity 
Class 6 cloud on a relative scale from 1 (least opaque) 
to 6 (most opaque). Lee & Myers (1999) identified 11 
cores within this complex, including L673-7, which they 
found to have a darkness contrast relative to the back- 
ground of 4 on a relative scale from 1 (lowest contrast) 
to 4 (highest contrast). They found that it was not as- 
sociated with an embedded Young Stellar Object (YSO) 
based on IRAS data and concluded that it was a starless 
core. Three related studies searching for infall motions 
towards starless cores using different dense gas tracers 
included L673-7 (Lee, Myers, & Tafalla 1999; Lee, My- 
ers, & Plume 2004; Sohn et al. 2007); none found any 
evidence for infall motions in this dense core. Park, Lee, 
& Myers (2004) included L673-7 in their survey of 12 CO 
J = 1— toward starless cores undertaken with the 14m 
Taeduk Ratio Astronomy Observatory (TRAO) in Ko- 
rea. They found that ^60% of their sample, including 
L673-7, showed some evidence for line wings that could 
result from outflows. They concluded that some of their 
starless cores may in fact harbor protostars, but higher 
angular resolution CO maps than the 50" achieved by 
their data and more sensitive infrared observations than 
available from IRAS were needed to determine which, if 
any, starless cores were in fact not starless. 

The distance to the L673 complex is somewhat uncer- 
tain. Herbig & Jones (1983) noted, based on optical im- 
ages, that it is clearly in the foreground of the extensive, 
large-scale Aquila obscuration, which begins at about 110 
pc and continues to about 1 kpc (Weaver 1949), setting 
an upper limit to the distance to L673. Herbig & Jones 
argue that the most likely distance to L673 is 300 pc 
based on astrometric methods. However, a careful read- 
ing of their discussion indicates that their results also 
allow for distances in the range of ~ 400 — 600 pc. In the 
following sections we adopt their most likely distance to 
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the L673 complex of 300 pc as the distance to L673-7. 
We discuss the effects of larger possible distances in Sj8] 

3. OBSERVATIONS 

3.1. Spitzer Space Telescope 

Observations of L673-7 were obtained with the Spitzer 
Space Telescope (Werner et al. 2004) in three different 
observing programs: Program ID (PID) 139 (c2d; Evans 
et al. 2003), PID 20386 (PI: P. C. Myers), and PID 30384 
(PI: T. L. Bourke). 

3.1.1. c2d 

Spitzer c2d (PID 139) observations of L673-7 were 
obtained with both the Infrared Array Camera (IRAC; 
Fazio et al. 2004) and the Multiband Imaging Photom- 
etry (MIPS; Rieke et al. 2004). The IRAC observations 
were obtained in all four bands, resulting in images at 
3.6 (IRAC band 1), 4.5 (IRAC band 2), 5.8 (IRAC band 
3), and 8.0 (IRAC band 4) fim. Two epochs of obser- 
vations were taken in order to identify and remove as- 
teroids. The first epoch was observed on 2004 October 
7 (AOR 0005152512); the second was observed on 2004 
October 8 (AOR 0005153024). Each epoch of observa- 
tions consists of two dithered, 12 s images. The dither 
size is approximately 10", resulting in an observed ficld- 
of-view approximately equal to the IRAC 5' x 5' field-of- 
view. The total integration time per pixel is 48 s. 

The MIPS observations were obtained in the first 
two photometry bands (24 and 70 /im). Two epochs 
were again observed to identify and remove asteroids. 
The first epoch was observed on 2004 April 14 (AOR 
0009427200); the second was observed on 2004 April 15 
(AOR 0009435904). The MIPS 1 (24 fim) data were 
obtained in the photometry mode using the small field 
size and a 3 s exposure time. A 1 column by 3 row raster 
map with full-array spacing in both the rows and columns 
was observed, giving an observed field-of-view of approx- 
imately 9'xl8' and a total integration time per pixel of 
84 s. The MIPS 2 (70 /im) data were also obtained in a 
nearly identical manner except 3 total cycles, rather than 
1, were observed per epoch, giving an observed field-of- 
view of approximately 4.5' x 13' and a total integration 
time per pixel of ~ 96 s. 

The IRAC and MIPS images were processed by the 
Spitzer Science Center (SSC), version S13, to produce 
basic calibrated data (BCD) images. The c2d Legacy 
Team then improved these images to correct for arti- 
facts remaining in the data, produced mosaics, performed 
source extraction, and compiled a band-merged source 
catalog. A complete description of the above steps can 
be found in Evans et al. (2007); see also Dunham et 
al. (2008) for a summary of the process. 

3.1.2. PID 20386 

Spitzer IRAC and MIPS data were also obtained in 
PID 20386 (PI: P. C. Myers), a program designed to ob- 
tain deeper observations of a sub-sample of the c2d cores. 
Hereafter we refer to these data as "cores2deeper" data. 
Two epochs of IRAC observations were obtained in all 
four bands. The first epoch was observed on 2005 Octo- 
ber 21 (AOR 0014606080); the second was observed on 
2005 October 24 (AOR 0014606336). Each epoch con- 
sists of 8 dithered, 30 s images, resulting in an approx- 
imately 5'x5' field-of-view with a total integration time 



per pixel of 480 s (10 times deeper in total integration 
time than the c2d observations). Additional 0.6 s "high- 
dynamic range" (HDR) mode images were obtained in 
each epoch to enable photometry on bright sources sat- 
urated in the longer exposures. 

The MIPS observations were again obtained at 24 and 
70 /im. One epoch was observed on 2005 September 27 
(AOR 0014614528). The 24 fim data were obtained in 
the photometry mode using the large field size with 3 cy- 
cles and a 10 s exposure time, giving an observed field-of- 
view of approximately 5.5'xl0.5' and a total integration 
time per pixel of 300 s (~ 3.5 times deeper in total in- 
tegration time than the c2d observations). The 70 /im 
data were obtained in an identical manner to the 24 /im 
data, giving an observed field-of-view of approximately 
3'xlO' and a total integration time per pixel of ~ 190 s 
(~ 2 times deeper in total integration time than the c2d 
observations) . 

The IRAC and MIPS images were processed by the 
SSC, version S14 for IRAC and S13 for MIPS, to produce 
BCD images. As with the c2d data, these cores2deeper 
BCD were then post-processed by the c2d pipeline (see 
references above). 

3.1.3. PID 30384 

Spitzer MIPS data were also obtained in PID 30384 
(PI: T. L. Bourke), a program designed to obtain MIPS 
band 3 (160 fim) images of select c2d cores. Hereafter 
we refer to these data as "coresl60" data. Due to the 
details of how MIPS operates, we gained additional 24 
and 70 /im data for free while obtaining the 160 /im data. 
The data were obtained in a single epoch on 2008 June 23 
(AOR 0018160896) in the scan map mode at the medium 
scan rate with 80" (~ 1/4 array spacing) return and for- 
ward leg cross scan steps. One cycle of 8 scan legs was 
obtained, resulting in a total integration time per pixel of 
160 s at 24 and 70 fim and 16 s at 160 fim. The 1/4 array 
spacing between scan legs ensured redundant coverage at 
160 /im. The data were reduced using the MIPS Data 
Analysis Tool (DAT; Gordon et al. 2005), as described 
in detail by Stutz et al. (2007). 

3.2. Caltech Submillimeter Observatory 

3.2.1. SHARC-II 350/450 fim Continuum 

Submillimeter continuum observations at 350 and 450 
/xm were obtained with the Submillimeter High Angular 
Resolution Camera II (SHARC-II; Dowell et al. 2003) 
at the Caltech Submillimeter Observatory (CSO). The 
350 fim observations were obtained in October 2008 
and will be fully described in an upcoming publication 
(M. M. Dunham et al. 2010, in preparation). The 450 
fim observations were obtained in June 2005 and are de- 
scribed in Wu et al. (2007). Details on the observing 
strategy, data reduction, and source photometry can be 
found in the above references. Table [TJ which lists wave- 
length, flux density, uncertainty in flux density, aperture 
diameter, and references for all available photometry on 
L673-7 (see below), presents flux densities at both wave- 
lengths. We list the 450 fim photometry for complete- 
ness; however, due to uncertainties introduced by uncer- 
tain calibration at this wavelength and loss of extended 
emission in the observing mode chosen for these obser- 
vations (Wu et al. 2007; M. M. Dunham et al. 2010, in 
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TABLE 1 
Photometry of L673-7 



A 


S V {X) 




Aperture 




H 


(mjy) 


(mjy) 


(arcsec) 


Reference a 


3.6 


<0.03 






1 


4.5 


<0.09 






1 


5.8 


<0.07 






1 


8.0 


<0.18 






1 


24 


1.30 


0.43 


6.0 b 


1 


70 


150 


26 


18 b 


1 


160 


<3800 






1 


350 


2700 


500 


40 


2 


450 c 


1100 


400 


40 


3 


1200 


480 


26 


120 


4 



a References - (1) New Spitzer observations; (2) 

M. M. Dunham ct al. (2010, in preparation) (3) Wu 

et al. (2007); (4) Kauffmann et al. (2008). 

b FWHM of Spitzer point-spread profile. 

c Not considered in all of the following analysis due to 

uncertain calibration and loss of extended emission. 

preparation), we omit this result from all analysis that 
follows. 

3.2.2. 12 CO J = 2—1 

Observations of 12 CO J = 2-1 at 230.537970 GHz 
were obtained in June and September 2006 at the CSO 
with the 230 GHz sidecab receiver and 50 MHz AOS 
backend (Kooi et al. 1998), providing 1024 channels with 
a resolution of 0.049 MHz per channel. The average in- 
strument resolution is approximately 2 channels per res- 
olution element, giving an average velocity resolution of 
0.13 km s -1 . The standard chopper wheel calibration 
method (Penzias & Burrus 1973) was used to measure 
T* A . The beam FWHM is 32.5" at 230 GHz 8 . The zenith 
optical depth at 225 GHz (measured via skydips every 
10 minutes) varied between approximately 0.05 — 0.07 
in June 2006 and 0.15 - 0.18 in September 2006, and 
the system temperature varied between approximately 
340 - 400 in June 2006 and 300 - 400 in September 2006. 

The main beam efficiency, ?/ m b, was calculated for the 
June 2006 data from observations of Jupiter. No plan- 
ets were available in September 2006, so ?y m b for this 
run was instead calculated by comparing observations of 
several bright sources to observations of the same sources 
obtained on previous runs with known efficiencies. The 
calculated f/ m b shows no systematic change between June 
and September 2006, with an average and standard de- 
viation of 0.77 and 0.08, respectively. All scans were 
converted from T\ to T m b (mean beam temperature) by 
dividing by f/ m b, and all temperatures given in this pa- 
per are in units of T m b. The pointing was checked and 
updated every 1 — 2 hours; the standard deviations in 
pointing were 6" in azimuth and 5" in zenith angle, giv- 
ing an overall pointing uncertainty of ~ 8". Since most 
of the pointing changes between updates were due to a 
slow drift with time, 8" is actually an upper limit to the 
true pointing uncertainty for any given observation. 

A nine-point, position-switched map with a spacing 
of 30" between scans (approximately full-beam spacing) 
was observed in June 2006. The map is centered at the 



See http://www.submm.caltech.edu/cso/receivers/beams.html 



position 19' l 21 m 34.9 s +11°21'19.0". A clean off-position 
of 8500" east (19 A 31 m 01.6 s +11°21'19.0") was found us- 
ing the 12 CO J = 1—0 Galactic Plane Survey (Dame 
et al. 2001) available through the Sky View Virtual Ob- 
servatory 9 . The map was then extended in June and 
September 2006, always with 30" spacing between scans, 
following the location of outflow emission seen in individ- 
ual scans (see Sj4T4} . A total of 86 distinct positions were 
observed, giving a total observed area of 21.5 arcmin 2 . 
The average la rms in the 86 scans is 0.25 K in the 0.13 
km s _1 channels. 

3.3. Other 

Millimeter continuum observations at 1.2 mm of L673- 
7 were obtained with the Max-Planck Millimetre Bolome- 
ter array (MAMBO) at the IRAM 30-m telescope. De- 
tails on the observing strategy and data reduction can be 
found in Kauffmann et al. (2008), and Table [1] lists the 
1.2 mm flux density of L673-7 taken from the photometry 
presented by Kauffmann et al. (2008). 

4. RESULTS 
4.1. Dark Core and Infrared Point Source 

Figures [T] and [2] show the Spitzer images of L673-7. 
Figure Q] shows two three-color images of L673-7: one 
comprised of IRAC 1 (3.6 /im; blue), IRAC 2 (4.5 /im; 
green) , and IRAC 4 (8.0 /im; red) c2d images, and one 
comprised of IRAC 4 (8.0 /mi; blue), MIPS 1 (24 /im; 
green), and MIPS 2 (70 /im; red) c2d images. Fig- 
ure [2] presents individually all of the IRAC and MIPS 
cores2deeper observations, along with SHARC-II 350 /im 
and MAMBO 1.2 mm continuum emission contours. 

Figures Q] and [2] show that the L673-7 core shows up in 
absorption against the mid-infrared background at 5.8, 
8.0, and 24 /im, as confirmed by the 350 /im and 1.2 
mm continuum emission contours tracing this dark ab- 
sorption region. Similar Spitzer dark cores have been 
presented for other low- mass dense cores (e.g., Stutz et 
al. 2007; Stutz et al. 2009a; Stutz et al. 2009b) and have 
often been called "shadows" ; here we instead adopt the 
term "dark core." 

An infrared point source, SSTc2d J192134.8+112123 
(J2000 HHMMSS.s+DDMMSS; hereafter referred to as 
L673-7-IRS) is visible at 24 and 70 /im in Figures [1] and 
[21 L673-7-IRS is not detected at any IRAC wavelength 
in the c2d band- merged source catalog except IRAC 2. It 
does have faint IRAC 1 — 3 detections (signal-to- noise ra- 
tios (S/N) between ~ 5 — 10) in the cores2deeper catalog. 
However, close visual inspection shows no clear evidence 
for a point source in the IRAC 1 — 3 images. There 
does appear to be very faint, extended, nebulous emis- 
sion in these images (not visible in Figures [1] and [5] due 
to the adopted image scalings) that may have confused 
the automated point-source extraction routine, but, if 
present at all, this extended emission is barely above the 
noise limit and likely arises from a combination of scat- 
tered light off the edges of an outflow cavity and emission 
from molecules excited by the molecular outflow ( ^4.41 
Noriega-Crespo et al. 2004). For the purposes of assem- 
bling a complete spectral energy distribution (SED) to 

9 Available at http://skyview.gsfc.nasa.gov/ 
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Fig. 1. — Left: Throe-color image of L673-7 comprised of the IRAC 1 (3.6 /im; blue), IRAC 2 (4.5 /im; green), and IRAC 4 (8.0 /im; 
red) c2d images. The color scales are displayed using a linear stretch in each band, with the minimum and maximum (0.3 and 2.7 MJy 
sr -1 for IRAC 1, 0.15 and 1.8 MJy sr -1 for IRAC 2, and 8.75 and 11.5 MJy sr^ 1 for IRAC 4, respectively) chosen to emphasize the dark 
core seen in absorption against the diffuse background at 8 /im. The white crosshair marks the position of L673-7-IRS. Right: Three-color 
image of L673-7 comprised of the IRAC 4 (8.0 /im; blue), MIPS 1 (24 fim; green), and MIPS 2 (70 fim; red) c2d images. The color scales 
are displayed using a linear stretch in each band, with the minimum and maximum (8 and 12 MJy sr — 1 for IRAC 4, 36 and 38 MJy sr 
for MIPS 1, and 60 and 100 MJy sr -1 for MIPS 2, respectively) chosen to emphasize both the point source at 24 and 70 fim and the dark 
core seen in absorption against the 8 and 24 fim diffuse backgrounds. The white crosshair marks the position of L673-7-IRS. Note that the 
edge of the 70 fim coverage is just to the west of L673-7-IRS, resulting in the black region of no data. 



calculate evolutionary indicators (see below) and con- 
strain radiative transfer models (see Sj6|, we treat the 
IRAC 1 — 3 cores2deeper detections as upper limits and 
list them as such in Table [TJ There is no IRAC 4 de- 
tection in either the c2d or cores2deeper data and no 
evidence for emission at the position of L673-7-IRS in 
the IRAC 4 images. We calculate a 3<r point-source sen- 
sitivity of 0.18 mJy for the IRAC 4 cores2deeper image, 
based on the standard deviation of the sky background, 
and list this as an upper limit in Table [T] 

L673-7-IRS is detected at 24 /xm in the c2d data and 
listed in the catalog with a flux density of >S^ 4Aim = 
1.94 ± 0.27 mJy (where the uncertaintly includes a 4% 
calibration uncertainty [Evans et al. 2007]). There is no 
24 /im detection listed in the cores2deeper catalog despite 
the clear existence of a point source in Figured We thus 
performed aperture photometry to calculate the 24 /im 
cores2deeper flux density of L673-7. We adopted an aper- 
ture radius of 7", sky annuli of 7"and 13", and an aper- 
ture correction of 2.05 as listed on the Spitzer Science 
Center website 10 . To properly calibrate our aperture 
photometry with the c2d and cores2deeper catalog PSF 
photometry, we performed aperture photometry with the 
same aperture radius, sky annuli, and aperture correction 
as above on the 50 sources in the cores2deeper catalog 
with a 24 /im detection with S/N > 3 and determined 
an additional correction factor of 0.88 by comparing the 
results to the catalog flux densities. With the above pa- 
rameters, we calculate S^ 4Aim — 1.32 ± 0.42 mJy (again 
including a 4% calibration uncertainty). 

L673-7-IRS is also detected at 24 /im in the coresl60 
data. We performed aperture photometry using the same 
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method as described above, including an additional cor- 
rection factor of 0.82 based on comparison of aperture 
photometry on the coresl60 data to cores2deeper cata- 
log photometry. We calculate S^" 1 = 1.10 ± 0.15 mJy 
(where the uncertainty again includes a 4% calibration 
uncertainty). We take a weighted average of the above 
three results and obtain a final S^ 4pm = 1.30 ±0.43 mJy, 
where the uncertainty is the standard deviation of the 
three separate measurements. We list this result in Ta- 
bled! 

At 70 /xm, the c2d catalog lists a flux density of 
gJOfim _ -j- 35 mJy (where the uncertaintly includes 
a 20% calibration uncertainty [Evans et al. 2007]). Auto- 
matic 70 /im source extraction was not performed on the 
cores2deeper data for technical reasons, and no source 
extraction was performed on the coresl60 data. Per- 
forming aperture photometry using a 16" radius aper- 
ture, sky annuli of 18" and 39", and an aperture cor- 
rection of 2.07 (as listed on the SSC website), we cal- 
culate Sl°^ m = 129 ± 30 mJy for the cores2deeper data 
and Sl°» m = 181 ± 37 mJy for the coresl60 data (again 
including 20% calibration uncertainties). No additional 
corrections are included because there are no other bright 
70 fim sources in the field for which such corrections can 
be obtained. A weighted average of the above values 
gives a final S™^ m = 150 ± 26 mJy, as listed in Table [U 

Figure H plots S^ m and Sl°^ m versus Julian Date for 
the multiple epochs of observations to search for evidence 
of variability. There is some evidence for a systematic de- 
crease in Sl^™ versus time over the four years spanning 
from April 2004 through June 2008. However, given the 
sparse sampling (only four epochs, two of which are sep- 
arated by only one day) and uncertainties in additional 
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Fig. 2. — Spitzer cores2deeper images of L673-7. Top Left: 3.6 (im image displayed in a linear stretch with the minimum and maximum 
intensities set to —3 and 5 Mjy sr — , respectively. The white crosshair marks the position of L673-7-IRS. Top Middle: 4.5 /an image 
displayed in a linear stretch with the minimum and maximum intensities set to —3 and 5 MJy sr — 1 , respectively. The white crosshair 
marks the position of L673-7-IRS. Top Right: 5.8 (im image displayed in a linear stretch with the minimum and maximum intensities set 
to 2 and 4 MJy sr — x , respectively. The white crosshair marks the position of L673-7-IRS. Bottom Left: 8.0 fim image displayed in a linear 
stretch with the minimum and maximum intensities set to 8.5 and 13 MJy sr -1 , respectively. The white crosshair marks the position of 
L673-7-IRS. The white contours show the SHARC-II 350 (im continuum emission. The contour levels are [4, 6, 10]rr, where la = 65 mjy 
beam -1 . The beamsize of the 350 fim observations is shown in the bottom right of the panel. Bottom Middle: 24 (im image displayed 
in a linear stretch with the minimum and maximum intensities set to 25.25 and 27 MJy sr -1 , respectively. The white contours show the 
MAMBO 1.2 mm continuum emission. The contour levels are [5, 7, 10, 13, 16, 19]rr, where la = 1.3 mjy beam -1 . The beamsize of the 1.2 
mm observations is shown in the bottom right of the panel. Bottom Right: 70 fim image displayed in a linear stretch with the minimum 
and maximum intensities set to 50 and 71 MJy sr — 1 , respectively. The contours are the same as in the previous panel. 

photometry correction factors derived by comparing cat- 
alog and aperture photometry flux densities, the case for 
this decrease remains ambiguous. There is no evidence 
for variability at 70 /xm within the uncertainties. 

Figure 0] shows the 160 /im image of L673-7 obtained 
in the coresl60 program. There is no clear detection 
associated with L673-7-IRS. Using the IDL Astronomy 
User's Library 11 program sky.pro, which calculates the 
sky background and variation in sky background using 
a procedure that iteratively eliminates outlier pixels, we 
calculate a sky background and variation in this back- 
ground of -8.0 MJy sr -1 and 23.3 MJy sr" 1 , respectively. 

Approximating the 160 /im Spitzer PSF with a Gaus- 
sian with #fwhm equal to the diffraction limited resolu- 
tion of ~ 40", this background intensity is then converted 
into units of mjy beam _1 as follows: 
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This gives the la point source sensitivity since the total 
flux density of a point source is equivalent to its flux 
density per beam. Substituting in the above value of 
er sky = 23.3 MJy sr -1 and multiplying by 3 gives a 3a 
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Fig. 3.— S^ Mm (top) and 5^ U/Jm (bottom) vs. time for the 
multiple epochs of MIPS observations spanning from April 2004 
through June 2008. 



11 Available at |http: / /idlastro.gsfc.nasa.gov/ 



point source sensitivity of 3000 mjy, which we list in 
Table Q] We consider this to be a conservative upper 
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Fig. 4. — Spitzer 160 /im coresl60 image, displayed in a linear 
stretch with the minimum and maximum intensities set to —75 
and +75 MJy sr _1 , respectively. The white crosshair marks the 
position of L673-7-IRS. 

limit since it encompasses large-scale sky variation over 
the map on top of the true local variation at the position 
of L673-7. 

Before moving on, we note here that there are 7 sources 
in the c2d source catalog for L673-7 classified as candi- 
date young stellar objects based on the selection criteria 
developed by Harvey et al. (2007). L673-7-IRS is not 
among these 7 sources because it is not detected at all 
Spitzer wavelengths between 3.6 — 24 /im, a requirement 
imposed by the Harvey et al. selection criteria. These 
7 sources have 2 — 24 /im spectral indices (see Evans et 
al. [2007] for details on the calculation of these indices) 
ranging from —2.52 to —1.82 and are thus classified as 
Class III objects (Greene et al. 1994; Evans et al. 2009). 
None are clearly associated with L673-7 itself; they are 
scattered about the ~ 5' x 5' field-of-view, with the clos- 
est of the 7 located 100" from L673-7-IRS (a projected 
separation of 30,000 AU at 300 pc). All 7 sources are 
classified as AGB stars according to the color and mag- 
nitude criteria developed by Robitaille et al. (2008) to 
perform a rough separation of YSO and AGB stars. Fi- 
nally, L673-7 is located close to the galactic plane in the 
direction of the inner galaxy (I ~ 46.5°, b ~ —1.5°), and 
Robitaille et al. (2008) found that the surface density of 
AGB stars increases exponentially with both decreasing 
galactic longitude and decreasing galactic latitude. Thus, 
we conclude that these 7 sources classified as candidate 
YSOs are in fact likely background AGB stars. 

4.2. Evolutionary Status of L673-7-IRS 

L673-7-IRS features a rising SED in the infrared (Table 
[T]), it is located within a dark absorption core in the mid- 
infrared and associated with submillimeter and millime- 
ter continuum emission, and it drives a molecular outflow 
(see N4.4p . Thus, L673-7-IRS is an embedded protostar, 
based on the criteria given by Dunham et al. (2008). Fig- 
ure[S]plots the observed SED of L673-7, including Spitzer 
MIPS detections at 24 and 70 ^m, the CSO SHARC-II 
detections at 350 /im, and a MAMBO 1.2 mm detection. 
Also shown are Spitzer IRAC 3.6-8.0 /tm and MIPS 160 
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Fig. 5. — Observed SED of L673-7. Detections are plotted as 
filled circles with error bars and upper limits are plotted as filled 
inverted triangles. See the text for details on photometry at each 
wavelength. 

/im upper limits . Based on the 24, 70, 350, and 1200 /im 
detections we calculate L^oi — 0.18 L Q , T^oi = 16 K, and 
Lboi/ Lsmm = 7 (where L smm is defined to be the lumi- 
nosity longward of 350 /im). The integrals are evaluated 
using the trapezoid rule (see Dunham et al. [2008] for de- 
tails and a discussion of errors introduced by integrating 
over finite, incompletely sampled SEDs). Based on both 
the calculated Tboi and L^oi/ L smrn , L673-7 is classified 
as a Class object (Andre et al. 1993; Chen et al. 1995). 

The calculated L^oi is likely a lower limit to the true 
bolometric luminosity since the SED is not fully sam- 
pled near its peak and the photometry at A > 100 /tm 
is performed in apertures likely to be smaller than the 
total source extent. Furthermore, the calculated L^oi is 
not a good indicator of the total Li nt since most of the 
bolometric luminosity likely arises from external heating 
by the interstellar radiation field (ISRF), which can add 
up to a few tenths of a solar luminosity to the bolomet- 
ric luminosity of an embedded, low- mass protostar (e.g., 
Evans et al. 2001). We can obtain an estimate of Li nt 
from the relation between Li nt and 70 //m flux found by 
Dunham et al. (2008). Applying their relation gives Li nt 
= 0.04 L . This estimate is good only to within about a 
factor of two. Radiative transfer models are required to 
better constrain both L^oi and Li nt and are presented in 

m 

4.3. Core Mass 

An estimate of the mass of the L673-7 core can be 
obtained from the 1.2 mm MAMBO data. Assuming 
that the 1.2 mm continuum emission is optically thin 
and the core is isothermal, the total mass of the core is 
given as 

M = 100 ff" , (2) 

where S v is the total flux density, B v (To) is the Planck 
function at the isothermal dust temperature To, k v is 
the dust opacity, d = 300 pc, and the factor of 100 is the 
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Fig. 6. — CO J = 2—1 spectra at three positions observed 
within the initial nine-point map: (0", 0"), (+30", +30"), and 
(-30", -30"). The spectrum at (+30", +30") is offset by +3 
K. It clearly shows a wing of blue-shifted emission relative to the 
spectrum at (0", 0"). The spectrum at (-30", -30") is offset by 
—3 K. It clearly shows a wing of red-shifted emission relative to 
the spectrum at (0", 0"). The dotted horizontal lines indicate the 
zero-levels of each spectrum, and the dotted vertical lines mark ± 
3 km s relative to the core systemic velocity of 7.3 km s . 

assumed gas-to-dust ratio. For consistency with later 
sections (see 3D and H]), we adopt K^ 2min = 1.02 cm 2 
gin -1 , the dust opacity of Ossenkopf & Henning (1994) 
appropriate for thin ice mantles after 10 5 yr of coagula- 
tion at a gas density of 10 6 cm -3 (OH5 dust). With the 
flux density listed in Tabled] and an assumed Tjj = 10.5 
K (the derived kinetic temperature, T^ n , of the core from 
NH 3 observations; T. L. Bourke et al. (2010), in prepa- 
ration), Equation [2] gives an estimated core mass of 1.9 
M Q . We list this result in Tabled 

4.4. Molecular Outflow 

As stated in ij3.2.2[ we initially observed a nine-point 
map with 30" spacing centered on L673-7-IRS 12 . Evi- 
dence for a molecular outflow was seen in this initial 
map in the form of a wing of blueshifted emission to the 
northeast and a wing of redshifted emission to the south- 
west. This evidence is shown in Figure [6l which shows 
the 12 CO J = 2— 1 spectra at three positions within the 
initial nine-point map. We then extended the map to 
follow this apparent northeast-southwest outflow axis. 

Figure [7] shows channel maps of the CO data at a ve- 
locity resolution of 2 km s _1 . The core systemic velocity 
is taken to be vi sr — 7.3 km s _1 from observations of 
C 18 J = 2-1, N 2 H+ J = 3-2, and DCO+ J = 3-2 
emission (J. H. Chen et al. 2010, in preparation). The 
central three channels are dominated by ambient emis- 
sion, but an outflow is clearly detected in the higher ve- 

12 The position of L673-7-IRS and the center of the 12 CO J = 
2—1 map actually differ by 4.6" because the position used for the 
CO observations was derived from an early version of the c2d data 
processing pipeline. Since this offset is smaller than the point- 
ing uncertainty in the 12 CO J = 2—1 map (8") and significantly 
smaller than both the spacing between scans and the beamsizc (30" 
and 32.5", respectively), the slight position offset between the two 
is negligible. 
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Fig. 7.— 12 CO J = 2-1 channel maps of the L673-7 outflow with 
a velocity resolution of 2 km s~ x . The grayscale shows main-beam 
temperature in K, with the mapping between grayscale and value 
of T m |-, given in the legend on the right. The contours are drawn at 
2, 5, 10, 20, 35, and 50cr, where the average lc rms is 0.095 K in 
the 2 km s — 1 channels. The lowest two contours are drawn with a 
black line; all others are drawn in white. The cross in each panel 
marks the position of the protostar, and the dotted line encloses 
the total area mapped. 

locity channels. Most of the outflow emission is found 
at \v — w C ore| < 4 km s _1 (where v core is the core sys- 
temic velocity), with no outflow emission detected be- 
yond \v — v core \ ~ 8 km s _1 . 

Figure [8] shows integrated redshifted and blueshifted 
emission. There is no obvious outflow emission at the 
map center, but there is one prominent blueshifted peak 
to the northeast and a corresponding redshifted peak 
to the southwest. The total extent of each peak is 
^150" from the center, and the position angle of the 
axis connecting the blue- and redshifted emission is ap- 
proximately 55° (measured east [counter-clockwise] from 
north). The fact that this outflow shows a clear bipolar 
morphology with mostly, but not entirely, distinct blue- 
and redshifted lobes and no outflow emission at the cen- 
ter suggests a moderate, but slightly edge-on, inclination 
(i ~ 45 — 70° ; between cases 2 and 3 of Cabrit & Bertout 
[1986]). A more precise estimate of the inclination would 
require knowledge of the outflow cone opening angle, but 
extreme pole-on and edge-on inclinations can be ruled 
out (the former by the bipolar morphology and the lat- 
ter by the mostly distinct red and blue lobes). There are 
additional redshifted peaks along the same axis located 
farther to the southwest, extending to the edge of the 
map, but the presence of corresponding blueshifted peaks 
is unknown since the map does not extend far enough to 
the northeast. The full extent of the molecular outflow 
driven by L673-7-IRS thus remains unknown. 

5. DARK CORE: COLUMN DENSITY PROFILE AND 
ENCLOSED MASS 

Figures [T]and [2] clearly show a dark region in the 8 and 
24 /im images that corresponds with the 350 fira and 1.2 
mm dust continuum emission. These dark regions result 
from absorption of the mid-infrared background by the 
dense core L673-7 and can be used to study the column 
density profile and total enclosed mass of the core, as has 
been done for prestellar cores with both Infrared Space 
Observatory (ISO) data (e.g., Bacmann et al. 2000) and, 
more recently, Spitzer data (e.g., Stutz et al. 2007, 2009a, 
2009b). The following analysis is based on the method 
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TABLE 2 
L673-7 Core Mass Estimates 





Aperture Radius 


Mass 


Source 


(arcsec) 


(Mq) 


1.2 mm flux density 


60 


1.9 


8 fim dark core 


70 


0.6 - 1.9 


24 jim dark core 


70 


0.2 - 4.5 


Dust Model 


83.3 


2.25 
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Fig. 8. — 12 CO J = 2—1 integrated intensity contours of L673- 
7 showing blueshiftcd and red-shifted emission, overlaid on a 24 
(im greyscale image centered on L673-7-IRS. The blue contours 
are integrated from -2.0 to 4.0 km s — while the red contours are 
integrated from 10.0 to 16.0 km s — 1 . The contours start at 0.5 K 
km s _1 and increase by 0.5 K km s _1 . The dotted line encloses 
the total area mapped. 

of Stutz et al., and uses the cores2deeper data. 
The observed intensity at a given frequency, I v , is given 

by 

I v = lie— + If" , (3) 

where 1° is the background intensity, l[ 9 is the fore- 
ground intensity, and r„ is the total optical depth. 
Knowledge of I? 9 is important for deriving an accurate 
intensity profile. Underestimating l[ 9 will cause one to 
underestimate the optical depth and corresponding col- 
umn density, while overestimating l[ 9 will cause one to 
overestimate the optical depth and corresponding col- 
umn density. Stutz et al. (2007) estimated l[ 9 from the 
darkest pixels in their map. However, the darkest pixels 
in both the 8 and 24 [im maps of L673-7 are found at 
the center of the dark core. Thus, since we have no di- 
rect way of measuring l£ 9 , we calculate lower and upper 
limits to the true absorption profiles by assuming either 
that there is no foreground emission (l£ 9 = MJy sr- 1 ) 
or that all of the emission in the darkest pixel is fore- 
ground emission (i.e., that the background emission is 
completely absorbed in the darkest pixel). In the latter 
case the assumed l£ 9 is subtracted from the image before 
the analysis below is conducted. 
Figured shows the 8 /im intensity profiles under the 



two above assumptions, and Figure H0k shows the 24 
/im intensity profiles. The intensities were obtained by 
averaging over concentric annuli, each with radii of 2.5" 
at 8 /im and 5" at 24 /im, centered at the position of 
L673-7-IRS. The positions plotted on the x-axis show 
the mid-points of the annuli. This analysis assumes that 
the core is spherically symmetric. The error bars show 
the standard deviation of the mean for each annulus; 
note that the actual standard deviations are much larger 
since the core is not perfectly spherically symmetric (see 
Figures [Hand 0). 

Inspection of the 8 and 24 /im intensity profiles shows 
that both profiles flatten out at ~70", are relatively flat 
between ~ 70 — 100", and then decrease again beyond 
100". Thus we define the edge of the core as traced by 
the absorption profile to be at a projected distance of 70" 
from the core center (corresponding to 21,000 AU at a 
distance of 300 pc), and average the values of J„ between 
70 — 100" to obtain the background intensities 1°. 

With 1° in hand, the optical depth through each an- 
nulus, t v , is calculated as 




(4) 



The 8 and 24 /im optical depth profiles are shown in 
Figures IHb andllOb. respectively. From the definition of 
optical depth as t„ = J n v p{r)dr, the total gas column 
density in each annulus is then 

N = 100 — — — , (5) 

where k v is the dust opacity, \i is the mean molecular 
weight per free particle (jjl = 2.37 for a gas that is 71% by 
mass hydrogen, 27% helium, and 2% metals; Kauffmann 
et al. 2008), toh is the mass of hydrogen, and the factor 
of 100 is the assumed gas-to-dust ratio. 

For consistency with the dust radiative transfer models 
in fjni we adopt the dust opacities of Ossenkopf & Hen- 
ning (1994) appropriate for thin ice mantles after 10 5 
yr of coagulation at a gas density of 10 6 cm~ 3 (OH5 
dust). The OH5 dust opacities give = 1063.8 cm 2 

g _1 and K 24Mm = 852.2 cm 2 g" 1 . For comparison, the 
Weingartner & Draine (2001; hereafter WD01) R v = 5.5 
dust opacities are about a factor of 1.5 lower than the 
OH5 opacities at 8 and 24 /im. However, Chapman et 
al. (2009a, 2009b) have shown that the mid-infrared ex- 
tinction law in dense regions flattens considerably more 
than accounted for by the WD01 Ry = 5.5 dust model 
and that the actual opacities may be up to a factor of 2 
higher than the WD01 Ry = 5.5 opacities at these wave- 
lengths. Thus we adopt the higher OH5 dust opacities 
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Fig. 9. — Analysis of the 8 fim dark core. In all three panels the 
filled circles with error bars show the results from assuming no fore- 
ground emission [jl a = MJy sr" 1 ) while the open squares with 
error bars show the results from assuming that all of the emission 
in the darkest pixel is foreground emission, (a): Intensity pro- 
file, calculated by averaging over concentric annuli, each with radii 
of 2.5", centered at the position of L673-7-IRS. The dashed and 
dotted lines show the background intensities obtained as described 
in the text for the assumptions of zero and maximum foreground 
intensity, respectively, (b): Optical depth profile, calculated from 
Equation [11 {left axis), and column density profile, calculated from 
Equation [5] {right axis). The solid and dotted lines show the col- 
umn density profiles for the best-fit p = 1.5 and p = 2.0 models, 
respectively (see JS}. (c): Total enclosed mass, calculated from 
Equation \E\ The solid and dotted lines show the total enclosed 
mass for the best-fit p = 1.5 and p = 2.0 models, respectively (see 
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Fig. 10. — Same as Figure [9] except now for the 24 jim dark 
core. In all three panels the filled circles with error bars show the 
results from assuming no foreground emission 

(j/s = o MJy sr- 1 ) 

while the open squares with error bars show the results from as- 
suming that all of the emission in the darkest pixel is foreground 
emission, (a): Intensity profile, calculated by averaging over con- 
centric annuli, each with radii of 5", centered at the position of 



L673-7-IRS. The left y-axis scale corresponds to the I, 



la 



MJy 



sr 1 intensity profile while the right y-axis scale corresponds to the 
il 3 = 24.9 MJy sr -1 intensity profile. The dashed and dotted lines 
show the background intensities obtained as described in the text 
for the assumptions of zero and maximum foreground intensity, re- 
spectively, (b): Optical depth profile, calculated from Equation [4] 
{left axis), and column density profile, calculated from Equation [5] 
{right axis). The solid and dotted lines show the column density 
profiles for the best-fit p = 1.5 and p = 2.0 models, respectively 
(see j|BJ. (c): Total enclosed mass, calculated from Equation [5] 
The solid and dotted lines show the total enclosed mass for the 
best-fit p = 1.5 and p = 2.0 models, respectively (see £|6)l . 
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but acknowledge there is about a factor of 2 uncertainty 
in these values. With the adopted dust opacities the to- 
tal gas column density in each annulus is calculated from 
Equation [5] and the results from the 8 and 24 /jm dark 
cores are shown in Figures [9Jd and 1 10b . respectively. 
The mass in each annulus is then calculated as 

M = (imnNA , (6) 

where A is the physical area of the annulus calculated 
assuming a distance of 300 pc. Figures andfTUb show 
the total enclosed mass calculated from the 8 and 24 fim 
dark cores, respectively. Based on the 8 /im dark core, 
the lower and upper limits to the total core mass are 0.6 
and 1.9 M©. Based on the 24 /im dark core, these limits 
are 0.2 and 4.5 M . These results are summarized in 
Table [2] 

Also plotted on Figures |9] and [10] are the column den- 
sity and total enclosed mass profiles from the best-fit 
p = 1.5 and p = 2.0 dust radiative transfer models. A de- 
scription of these models, along with a discussion of how 
well they match the column density and enclosed mass 
profiles calculated from the 8 and 24 /zm dark cores, is 
presented in Sj6] 

6. DUST RADIATIVE TRANSFER MODELS 

In order to obtain an accurate constraint on the inter- 
nal luminosity of L673-7-IRS, we have constructed phys- 
ical models of this source, constrained primarily by the 
source SED (Figure [5]) but also by the 350 /im and 1.2 
mm intensity profiles. The modeling procedure follows 
closely that presented by Dunham et al. (2006) (see also 
Young et al. 2004; Bourke et al. 2006; Lee et al. 2009; 
Young et al. 2009; Lee et al. 2010). 

We used the two-dimensional, axisymmetric, Monte 
Carlo dust radiative transfer code RADMC (Dullemond 
& Turolla 2000; Dullemond & Dominik 2004) to calcu- 
late the temperature profile and emergent SED of a pro- 
tostar embedded in a dense core (the term envelope is 
often used interchangeably with dense core in the lit- 
erature and in this paper). Despite RADMC allowing 
for two-dimensional structure, we only consider spheri- 
cally symmetric models to avoid introducing too many 
free parameters into models constrained by only a lim- 
ited number of data points (as given in Tabled]). For the 
dust properties, we adopt the dust opacities of Ossenkopf 
& Henning (1994) appropriate for thin ice mantles af- 
ter 10 5 yr of coagulation at a gas density of 10 6 cm~ 3 
(OH5 dust), which have been shown to be appropriate 
for cold, dense cores (e.g., Evans et al. 2001; Shirley et 
al. 2005). Isotropic scattering off dust grains is included 
in the model as described by Dunham et al. (2010). 

The envelope is heated both internally by the proto- 
star and externally by the ISRF. The external heating 
is included in our model; we adopt the Black (1994) 
ISRF, modified in the ultraviolet to reproduce the Drainc 
(1978) ISRF. Most of the models of low-luminosity, 
embedded protostars presented in the references listed 
above then attenuate this Black-Draine ISRF by some 
Ay of dust with properties given by Draine & Lee (1984) 
in order to provide the correct amount of external heating 
to match the peak of the observed SED. The physical jus- 
tification for this attenuation is to simulate the core being 
embedded in a parent cloud. For our models of L673-7 



any such attenuation causes the model to under-predict 
the 350 /zm emission. Increasing the internal luminos- 
ity (see below for more details on the internal source) 
enables the model to match the observed 350 /im emis- 
sion, but over-predicts the amount of 3.6 — 70 /im emis- 
sion. The 3.6 — 24 /zm model emission can be adjusted by 
changing the envelope inner radius (again, see below for 
more details), but the 70 /im emission, which primarily 
traces internal luminosity (Dunham et al. 2008), remains 
over-predicted. Thus, we do not attenuate the ISRF. Ul- 
timately there is a degeneracy between attenuation by 
surrounding material and intrinsic strength of the ISRF 
that cannot be broken by SED modeling, particularly an 
SED as sparsely sampled near the peak as our SED of 
L673-7. Observations at 100 — 500 /im with, for exam- 
ple, the Herschel Space Observatory, will be required to 
better constrain the peak of the SED and total amount 
of external heating. 

With the ISRF constrained as described above, the 
mass of the envelope is constrained to be M env ~ 2.25 
M Q by requiring that the model flux match the optically 
thin 1.2 mm continuum flux. We assume a power-law 
radial density profile [n(r) oc r~ p ] for the envelope and, 
following Young et al. (2009), we run two model grids, 
one with p = 1.5 and one with p = 2.0. Estimates of the 
total core size range from —10,000 AU from the SHARC- 
II 350 /im data to -40,000 AU from the MAMBO 1.2 mm 
data, with the 8 and 24 /im absorption profiles suggesting 
an intermediate value of —21,000 AU. We adopt an en- 
velope outer radius of 25,000 AU because it is both the 
midpoint between the SHARC-II and MAMBO source 
sizes and close to the core size derived from the absorp- 
tion profiles, but note that the true core size remains 
uncertain. The envelope inner radius (t^ v ) is allowed to 
vary between 1 - 300 AU as 1, 5, 10, 20, 40, 60, 300 
AU, for a total of 18 different inner radii. 

The internal source consists of a protostar and a cir- 
cumstellar disk. The protostar is assumed to be a 3000 
K blackbody; as noted by Dunham et al. (2006), these 
models are not sensitive to the protostellar temperature 
since most of the protostellar emission is reprocessed to 
long wavelengths. The protostellar luminosity (L star ) is 
allowed to vary between 0.0005-0.04 L as 0.0005, 0.001, 
0.005, 0.01, 0.015, 0.04 L , for a total of 10 different 
Lstar- The presence of a disk is inferred from the presence 
of a molecular outflow ( §4.4| and is in fact required for 
the model to match the 24 /im flux. Following Dunham 
et al. (2006), we include a disk in our 1-D model using 
the method developed by Butner et al. (1994), based on 
the disk model of Adams et al. (1988). This method 
simulates a disk by calculating the emission from a disk 
with given surface density and temperature profiles at 
a given inclination, averaging the emission over all in- 
clinations, and then adding this average emission spec- 
trum to the protostellar spectrum to form the final input 
SED of the internal source used to calculate the radiative 
transfer through the envelope. Both the surface density 
and temperature profiles are described as power laws: 
X(R disk ) oc i?"™ and T(R dlsk ) cx R^ k . We set w = 
1.5, following Butner et al. (1994) and Chiang & Goldrc- 
ich (1997), and q — 0.5 to simulate a flared disk. The 
model also includes an intrinsic disk luminosity, Ldi S k, 
perhaps generated by accretion within or onto the disk. 
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Ldisk varies over the same range as above for L star . 

The inner radius of the disk is set to be the radius at 
which the temperature is equal to the dust destruction 
temperature, calculated assuming spherical, blackbody 
dust grains and a dust destruction temperature of 1500 
K. The outer radius is set to the centrifugal radius. Given 
the sound speed, c s , the age of the object, t, and the 
angular velocity of the cloud prior to collapse, flo, the 
centrifugal radius is calculated as Rc = (mQ7l6)c s t 3 Qo, 
where mo is a dimensionless constant of order unity 
(Terebey et al. 1984; Young & Evans 2005). Assuming a 
gas temperature of 10 K, a sound speed calculated from 
equal thermal and turbulent contributions, t — 6.1 x 10 4 
yr (outflow dynamical time; see 37. 1[) . and Hq — 4 x 10~ 13 
s" 1 (Andre et al. 1999), we calculate R c = 128 AU. For 
all models with r™^ > 128 AU, the disk outer radius 
is set to Rc] otherwise, the disk outer radius is set to 
r™ v . The disk mass is assumed to be Mdisk = 0.01 M . 
Neither the disk outer radius nor the disk mass are con- 
strained by these models since the long- wavelength emis- 
sion detected by the MAMBO data is dominated by the 
envelope. Sensitive submillimeter and/or millimeter con- 
tinuum interferometer observations are required to study 
the disk size and mass (e.g., J0rgensen et al. 2009; Enoch 
et al. 2009b). 

With the above assumptions, we ran a grid of models 
over the two possible values of p, 18 possible values of 
r™^, 10 possible values of L sta r, and 10 possible values 
of Ldisk, for a total of 3600 models, 1800 each for the p = 
1.5 and p = 2.0 density profile power-law indices. Each 
model is compared to the observed SED by calculating a 
reduced chi-squared (xr) as 

2 _ 1 » [S o bS{Xi) _ Srd{Xi)] 2 

where S° bs (Ai) is the observed flux density at \, 
S™ od {\i) is the flux density from the model at this wave- 
length (calculated in an aperture size identical to that 
used to measure the photometry), and a u (\i) is the un- 
certainty in the observed flux density. For n data points 
and m free parameters there are k = n — m degrees of 
freedom in these models, and we divide by k to obtain 
the reduced \ 2 values. There are four data points (ex- 
cluding the 450 /im point) and three free parameters for 
each of the p = 1.5 and p = 2.0 model grids (r^ nv , L sta r, 
and Ldisk), resulting in one degree of freedom (k = 1). 
Finally, we check whether each model is consistent with 
the 3.6 — 8.0 and 160 fim upper limits 13 ; only models 
consistent with these upper limits are kept for further 
analysis. 

Figure Qj] shows the results from the p = 1.5 model 
grid. Figures [TTa and lllb show the minimum value of Xr 
for each value of L int (L int = L star + Ldisk) in the model 

13 The 160 fim upper limit listed in TablelTlis an upper limit for a 
point source; an extended source could have a higher total flux and 
still be undetected. To compare the models with this upper limit, 
we calculate the model flux within an aperture with a radius equal 
to twice the standard deviation of the 160 /an PSF and compare 
this result to the derived upper limit, since nearly all (95%, to be 
precise) of the flux from a point-source would be contained within 
this aperture. Assuming the PSF is a Gaussian with a FWHM of 
40", this radius is 34". 




50 100 150 200 250 300 10 20 30 40 50 

Envelope Inner Radius (AU) Envelope Inner Radius (AU) 

Fig. 11. — Results from the p = 1.5 model grid, (a): Minimum 
value of Xr f° r each value of L;„t in the model grid, (b): Same 
as (a), except with different x and y axis ranges. The dotted line 
shows the 99.9% confidence limit of 10.827. (c): Minimum value 
of Xr f° r each value of r™„ in the model grid, (d): Same as (c), 
except with different x and y axis ranges. The dotted line shows 
the 99.9% confidence limit of 10.827. 

grid. We do not show results for L sta r and Ldisk sepa- 
rately because, as noted by Dunham et al. (2006), these 
models are sensitive to the total Li nt but not particularly 
sensitive to the individual stellar and disk contributions. 
From Figure ITTk it is clear that the best- fit models have 
Lint ~ 0.02 Lq, and FigurefTTb demonstrates that impos- 
ing a 99.9% confidence limit of xf- < 10.827 constrains 
0.01 < L int (L ) < 0.03. Figures Eb and [TTJi, which 
plot the minimum value of Xr for each value of r^ v in 
the model grid, show that a similar 99.9% confidence 
limit constrains r|™„ < 5 AU. This constraint on the in- 
ner envelope radius should not be over-interpreted; it is 
simply the r*™„ that produces the correct total amount of 
optical depth through the model to match the observed 
SED, given the assumed density profile. Different den- 
sity profiles would require different r" nv to give the same 
total opacity (see below). 

Figure [12] is similar to Figure [TTJ but now shows the 
results for the p = 2.0 model grid. The 99.9% confidence 
limits constrain 0.035 < L mt (Lq) < 0.045 and 100 < 
r™. u (AU) < 150. Again, the constraint on the inner en- 
velope radius should not be over-interpreted. The con- 
straints on the internal luminosity, however, are robust, 
as demonstrated by the similar constraints independent 
of assumed density profile and as discussed by Dunham 
et al. (2006). 

Figure [TJ shows the SED of L673-7 with the best-fit 
models from the p = 1.5 and p = 2.0 grids overplotted. 
The best-fit p = 1.5 model has Li nt = 0.02 L (L s tar = 
0.015 L , L dlsk = 0.005 L ), r*£„ = 5 AU, and xl = 
3.5, while the best-fit p = 2.0 model has Li nt = 0.04 
L (Lstar = 0.01 L , L dlsk = 0.03 L ), rf m = 140 
AU, and Xr = 8.3. Both models are consistent with the 
3.6 — 8.0 /Ltm and 160 /im upper limits. The p = 1.5 
model appears to provide a slightly better fit by eye, a 
fact that is confirmed by the lower Xr, and also appears 
more consistent with the upper limits (although, strictly 
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Fig. 12. — Same as Figure [TT] except now showing the results 
from the p = 2.0 model grid. 

speaking, no upper limits are violated by either model). 
In general, however, we are unable to distinguish between 
these two models based on the SEDs alone. 

Unlike the SEDs, the submillimeter and millimeter in- 
tensity profiles are sensitive to the envelope density pro- 
file (e.g., Shirley et al. 2002). Figure [HI plots the model 
350 and 1.2 mm intensity profiles for the best-fit 
p = 1.5 and p = 2.0 models shown in Figure [T31 along 
with the observed intensity profiles calculated from the 
data in a manner essentially identical to the calculation 
of the 8 and 24 /Ltm intensity profiles in fJS] The model 
intensity profiles for the p = 1.5 model are clearly too 
shallow compared to the data. On the other hand, the 
p = 2.0 model intensity profiles provide a much better 
match to the observed profiles. 

Finally, we overplot the column density and enclosed 
mass profiles calculated from the best- fit p = 1.5 and 
p = 2.0 models on Figures [9] and [TQl We refrain from 
making any quantitative comparisons between the model 
and absorption profiles since we were only able to set 
lower and upper limits to the absorption profiles. Quali- 
tatively, the shapes of the model profiles generally match 
those derived from the 8 and 24 /im absorption profiles. 
The model profile shapes do appear to diverge from the 
absorption profiles beyond ~ 40" (12000 AU at 300 pc) 
in that the model column density profiles fall off less 
steeply and the model enclosed mass profiles become less 
flat than the absorption profiles. This divergence can 
easily be explained by the fact that the models assume 
an outer envelope radius of 25,000 AU, while the ab- 
sorption profiles flatten out to the assumed background 
intensity at 21,000 AU. These same best-fit models mod- 
ified to have outer envelope radii of 21,000 AU do indeed 
provide better matches to the overall shapes of the ab- 
sorption profiles (not shown). 

Based on the above results, we conclude that L673-7 
features a relatively steep density profile, with p ~ 2. 
A more accurate constraint on the exact density profile, 
particularly at small radii, would require both better lim- 
its on the absorption profiles (through better knowledge 
of the foreground intensities) and quantitatively testing 



more than just two density profiles, and is beyond the 
scope of this paper. From the fits to the SED alone 
we find that 0.01 < L int (L ) < 0.045. Once the con- 
straint on the density profile is added in, we conclude 
that the most likely value is L;„ t ~ 0.04 L Q . This is 
in very good agreement with the result from Section [4.21 
that Li n t — 0.04 Lq, based on the Dunham ct al. (2008) 
relation between Li nt and 70 /im flux. 

7. OUTFLOW PROPERTIES AND PROTOSTELLAR MASS 
ACCRETION RATE 

7.1. Outflow Properties 

We calculate the mass of the outflow, Mfi ow , and 
its kinematic and dynamic properties (kinetic energy, 
Efiow, momentum, Pfi ow , luminosity, Lfi ow , and force 14 , 
Ffiow) m the standard manner (e.g., Cabrit & Bertout 
1990). To calculate the mass (and all the other prop- 
erties that depend on mass), we first calculate the col- 
umn density of H2 , Ah, , within each velocity channel in 
each pixel. Assuming LTE, A H2 = f{T,X co )(J T mb dv), 
where f(T, Xco) is a function of the temperature of out- 
flowing gas, T, and CO abundance relative to H 2 , Aco- 
The integral J T m b dv is over the velocity channel and is 
given by the temperature in that channel multipled by 
the channel width (0.13 km s _1 ). Following Hatchell et 
al. (2007), we assume a standard CO abundance relative 
to H 2 of Aco = 10~ 4 . This abundance is generally un- 
certain to within a factor of ~ 3 (e.g., Frerking et al. 1982; 
Lacy et al. 1994), directly leading to a corresponding fac- 
tor of ~ 3 uncertainty in Mfi ow and all other quantities 
that depend on Mfi ow . 

The correct temperature to assume is unknown. 
Hatchell et al. (2007) chose 50 K, citing evidence for 
warm gas in outflows (e.g., Hatchell et al. 1999). Here we 
take a different approach. The temperature is only one 
of several sources of uncertainty in deriving the outflow 
properties. Ultimately, strong lower limits to these prop- 
erties are more useful than highly uncertain values that 
attempt to "average out" the various sources of uncer- 
tainty. Thus, we choose the temperature from the range 
of 10 — 100 K (encompassing the likely range for gas in 
molecular outflows) that minimizes /, 17.6 K. Choosing 
100 K instead, the temperature that maximizes /, would 
increase the calculated column density by a factor of 2.5. 
With T = 17.6 K, 



A H2 = 5.144 x 10 18 I J T mh dv) cm~ 2 . (8) 

The mass within each velocity channel in each pixel is 
then calculated as M VtP i xe i = ^h 2 to hAh 2 A P i xe i, where 
mn is the mass of a hydrogen atom, /ih 2 is the mean 
molecular weight per hydrogen molecule (pa 2 = 2.8 for 
gas composed of 71% hydrogen, 27% helium, and 2% 
metals by mass [Kauffmann et al. 2008]), and A p i xe i is 
the area of each pixel. The total mass of the flow is then 
obtained by summing M v . p i xe i over all velocity channels 

14 The outflow force is also commonly referred to as the outflow 
momentum flux. It is calculated as a momentum over time and 
therefore is in fact a force. Since the usual definition of the flux of 
some quantity is quantity per area per time, rather than quantity 
per time, we refrain from calling Ff[ ow the momentum flux in this 
paper. 
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Fig. 13. — Best-fit models of L673-7. The solid lines show the total mode l fluxes, the filled circles show the photometry listed in Table 
[T] (except for the 450 fim data point, which is not plotted; see text in £|3.2.1l for details), the filled inverted triangles show the upper limits 
listed in Table [T] and the open squares show the model fluxes in the aperture sizes used for photometry at long wavelengths (A > 100 
(im) where the total emission is more extended than the aperture. (Left): Best-fit model from the p = 1.5 grid, with Li n t = 0.02 Lq 
(L 3 tar = 0.015 L , L disk = 0.005 L ), r™„ = 5 AU, and xl = 3.5. (Right): Best-fit model from the p = 2.0 grid, with L mt = 0.04 Lq 
(L 3tar = 0.01 L Q , L dlsk = 0.03 L Q ), r*£„ = 140 AU, and x 2 r = 8.3. 
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Fig. 14. — Intensity profiles for L673-7. The solid and dashed lines show the 350 /i m a nd 1.2 mm model intensity profiles, respectively, 
for the best-fit models from the p = 1.5 (left) and p = 2.0 (right) grids shows in Figure [T3l The filled circles and squares show the observed 
350 (im and 1.2 mm intensity profiles, respectively. 

and all pixels. Following Hatchell et al. (2007), we in- 
clude all velocity channels for which \v — v core \ =2 — 15 
km s -1 (where v core is the core systemic velocity). Since 
we do not know the total extent of the outflow, we restrict 
the summation to those pixels within 150" of L673-7-IRS 
(the approximate extent of the main blue and red lobes; 
see H4.4I) and consider the result a lower limit to the to- 
tal outflow mass. Following this procedure, we calculate 
Mf^ > 0.05 Mq and list this result in Tabled 

The outflow kinetic energy and momentum within each 
velocity channel in each pixel are calculated as E VtP i xe i 



where v is the velocity of each channel (with respect to 
the core systemic velocity) . The total kinetic energy and 
momentum of the outflow are then calculated by sum- 
ming over all pixels and velocity channels, again includ- 
ing only those channels for which \v — v core \ = 2—15 
km s _1 and pixels within 150" of L673-7-IRS. We calcu- 



late E 



flow 



> 6.1 x 10 



42 1 

cos 2 



ergs and Pfi ow > 0.13 



\M v<p ixei x v 2 and P ViPixe i = M VjP i xe i X v, respectively, 



Mq km s^ 1 , and list these results in Table [3] The incli- 
nation dependence in both quantities arises from the fact 
that we measure radial velocities rather than true space 
velocities and is left in the results at this stage. 
The outflow luminosity and force are calculated as 
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TABLE 3 

Outflow and Mass Accretion Properties 



Quantity Units Value 

Outflow Mass (M f i ow ) M Q > 0.05 

Outflow Kinetic Energy (E flow ) 10 42 ergs > 6.1 

Outflow Momentum (Pfi aw ) M © km s -1 > 0.13 — 

Outflow Luminosity {L flow ) 10" 4 L Q > 8.3 s< ^ 

Outflow Force (F flow ) 10 -6 M km s" 1 yr -1 > 2.1 — - 

Outflow Force a (F flow ) 10~ 6 M km s" 1 yr^ 1 > 4.5 f^- 

Average Protostellar Mass Accretion Rate «M acc )) 10~ 6 M Q yr _1 > 1.2 

Dynamical Time (r<j) 10 4 yr > 6.1 



sin i 

Protostellar Mass Accreted (M aC c) M > 0.07 

Accretion Luminosity 13 (L acc ) Lq > 0.36 L© 



a Calculated using the force per beam method (Fuller &; Ladd 2002; Hatched ct al. 2007). 
k Assuming constant mass accretion so that average rate is equal to current rate. 



Lfiow = E f i ow /T d and F f i ow = Pfi ow /r d , respectively, 
where r d is the dynamical time of the outflow. Assum- 
ing a total extent of at least 150" and a characteristic 
velocity of v = Pfiow/Mfiow = 3.5 km s -1 , we calculate 
r d > 6.1 x 10 4 fff yr, L flow > 8.3 x 10 L 01 and 

Ffi ow > 2.1 x 10" 6 ^i M km s" 1 yr" 1 , and list these 
results in Table El °°K gain, we leave the inclination de- 
pendence in the results at this stage. These lower limits 
on the mass, luminosity, and force of this outflow are 
generally consistent with the correlations between these 
quantities and the luminosity of the driving source pre- 
sented by Wu et al. (2004), after extrapolating by eye 
these correlations (and their disperstions of several or- 
ders of magnitude) to the very low luminosity of L673-7- 
IRS. However, since we are only able to determine lower 
limits to the outflow properties, it is also possible that 
the L673-7 outflow has a higher mass, luminosity, and/or 
force for the luminosity of its driving source compared to 
other outflows. 

The outflow force is essentially the average rate at 
which momentum has been injected into the outflow over 
the region encompassed by the calculation of the outflow 
momentum and dynamical time (in this case, 150", the 
extent of the primary red and blue lobes). However, the 
above calculation simply divides the total momentum by 
the total dynamical time, and does not take into account 
the fact that the total time that outflowing gas has spent 
in the outflow is a function of the velocity of the gas 
(at the same distance from the protostar, lower velocity 
gas has been in the outflow for longer than higher ve- 
locity gas). We thus also calculate the force using the 
outflow force per beam method (Fuller & Ladd 2002; 
Hatchell et al. 2007), which calculates the contributions 
to the total force separately for each velocity channel and 
then sums over all velocities. Again including all veloc- 



ity channels for which \v 
calculate Fri ow > 4.5 x 10 



2 - 15 km s"\ 



-6 sin i 



M km s 



yr" 



We 



list this result in Table [H The two calculations of the 
outflow force agree to within a factor of about 2, but we 
consider this second calculation to be more accurate and 
use it in the discussion below. 

We have not included an opacity correction of 1 _ T c ~r 
in the calculation of the outflow mass (and thus in all 
other quantities that depend on the outflow mass) to 



account for the fact that the 12 CO J = 2— 1 emission 
is likely optically thick. This factor is always greater 
than or equal to 1, and is typically about 3.5 (Cabrit 
& Bertout 1992; Bontemps et al. 1996). This opacity 
correction, combined with the discussion above about the 
LTE excitation temperature, implies that our calculated 
properties are strong lower limits. 

7.2. Protostellar Mass Accretion Rate 

Molecular outflows such as the one presented here are 
driven by the transfer of momentum from a jet/wind 
ejected by the protostellar system to the ambient medium 
(e.g., Bontemps ct al. 1996; Richer et al. 2000; Arce et 
al. 2007). Following Bontemps et al. 1996, conservation 
of momentum requires that 



Fflow — fcntM w V w 



(9) 



where M w is the mass- loss rate in the jet/wind, V w is the 
velocity of the jet/wind, and / en t is the entrainment effi- 
ciency between the jet/wind and the ambient gas. This 
expression can be re-written as 



flow 



= fe 



Mw 

Mac 



V m M a , 



(10) 



where M acc is the protostellar mass accretion rate. 

The rate at which mass is lost in the jet/wind is cor- 
related with the protostellar mass accretion rate (Bon- 
temps et al. 1996 and references therein). The term 
Mu,/M acc in Equation [TU] can thus be considered a con- 
stant, resulting in a direct proportionality between the 
protostellar mass accretion rate and the outflow force: 



A/a 



1 



feat M w V w 



F 



flow 



(11) 



We assume a typical jet/wind velocity of V w ~ 150 
km s _1 (Bontemps et al. 1996). Models of jet/wind for- 
mation predict, on average, M w /M acc ~ 0.1 (Shu et 
al. 1994; Pelletier & Pudritz et al. 1992; Wardle & Konigl 
1993; Bontemps et al. 1996). The entrainment efficiency 
is typically f cnt ~ 0.1—0.25 (Bontemps et al. 1996; Andre 
et al. 1999). Since we calculated a lower limit to Ffi ow 
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and thus will also calculate a lower limit to M acc , we 
adopt / ent = 0.25. With these values and our result for 
F f i ow , we calculate (M acc ) > 1.2 x KT 6 ^ M Q yr' 1 , 
close to the canonical value of ~ 2 x 10~ 6 Mq yr -1 for the 
inside-out collapse of a singular isothermal sphere (Shu, 
Adams, & Lizano 1987). The brackets indicate that this 
is the time-averaged mass accretion rate over the lifetime 
of the outflow, since Ffi ow is the average rate at which 
momentum has been injected. 

With knowledge about the age of the outflow and 
the above result for the time-averaged protostellar mass 
accretion rate, we can calculate the protostellar mass 
accreted. For the age of the outflow we use the dy- 
namical time, Td > 6.1 x 10 4 |^| yr. The protostellar 

mass accreted is then M acc = M acc r d . We calculate 
M acc > 0.07 ^— M Q and list this result in Table d 

We now have a lower limit to both the time-averaged 
protostellar mass accretion rate and the protostellar 
mass. With these quantities we can calculate the ac- 
cretion luminosity, 

^ GM acc M acc . . 

Lace — „ 7 

H 

where G is the gravitational constant and R is the proto- 
stellar radius (assumed to be 3 Rq 15 ). With our results 
from above, we find that L acc > 0.9 . L . Through 
this point we have left the inclination dependence in our 
results. Since we are calculating a lower limit to the ac- 
cretion luminosity, we evaluate the expression e s o '" 3 ^ at 
the inclination that minimizes this quantity. We restrict 
ourselves to the range 20° < i < 70° since extreme pole- 
on and edge-on inclinations can be ruled out ( ij4.4[) . The 
minimum value of 0.41 occurs at i = 20°, while the max- 
imum value occurs at z=70° and is larger by a factor of 
57. Inserting i = 20° into the above expression for L acci 
we find 

L acc >O.36L . (13) 

We list this result in Table [3] 

Our result for the accretion luminosity is a strong lower 

2 

limit since L acc oc M acc M acc oc M acc r d and both M acc 
and Td are lower limits. Specifically, M acc can only in- 
crease (typically by a factor of ~ 3.5) when an opacity 
correction is included to account for the fact that the 
12 CO J = 2—1 emission is likely optically thick, it can 
only increase (by up to a factor of 2.5) if the outflowing 
gas is not at the assumed temperature of 17.6 K, and it 
can only increase (by up to a factor of 57) if the incli- 
nation is larger than 20° (which it likely is; H4.4I) . Ad- 
ditionally, the true dynamical time is almost certainly 
larger by at least a factor of 2 since the map presented 
here does not cover the full extent of the outflow. Thus, 
we conclude with strong confidence that L acc > 0.36 Lq, 
with the actual value likely much higher. 

5 A larger accretion radius, as would be the case if material were 
accreting only onto a disk, would lead to a correspondingly smaller 
accretion luminosity. Indeed, most material accreting from the 
envelope likely accretes onto the disk (e.g., Terebey, Shu, & Cassen 
1984; Young & Evans 2005). However, the presence of the outflow 
argues in favor of accretion from th e di sk onto the protostar, and 
since i? pro tostar << Kdiski Equation ll2l is a good approximation of 
the accretion luminosity generated by this accretion. 



8. DISCUSSION 

With the detection of an infrared source driving a 
molecular outflow, L673-7 is clearly a protostellar core. 
However, as noted in three surveys for infall motions 
towards starless cores included L673-7 but did not find 
any evidence for infall in this core (Lee, Myers, & Tafalla 
1999; Lee, Myers, & Plume 2004; Sohn et al. 2007). The 
molecular line observations used to search for infall signa- 
tures in these three studies featured beam FWHM rang- 
ing from ~ 25 — 65", corresponding to physical sizes of 
7500 — 19500 AU at the assumed distance of 300 pc. In 
the Shu (1977) standard model, a wave of infall propages 
outward at the sound speed, leading to the definition of 
an infall radius r; n f = c s t, where c s is the sound speed 
and t the time since collapse began. Material within 
the infall radius is collapsing onto the central protostar 
whereas material outside the infall radius is still at rest. 
Assuming a purely thermal sound speed with 10 K gas, 
the infall radius reaches radii of 7500 — 19500 AU at 
- 190,000 - 488,000 yr after collapse begins. These 
times are a factor of 3 — 8 larger than the outflow dynam- 
ical time. If the outflow dynamical time represents the 
age of the protostar then the collapsing region is smaller 
than the size of the beam, possibly explaining the non- 
detection of infall signatures. In reality the derived out- 
flow dynamical time is only a lower limit since we have 
not mapped the full extent of the outflow, and material 
outside the infall radius might already be infalling rather 
than at rest (e.g., Fatuzzo, Adams, & Myers 2004). How- 
ever, both rotation and outflows can mask the kinematic 
signatures of inall (Lee, Myers, & Tafalla 1999 and ref- 
erences therein) , thus a nondetection of infall signatures 
does not necessarily imply a lack of infall. 

Based on both a strong correlation between the 70 
fim flux of a protostar and its internal luminosity (Dun- 
ham et al. 2008) and dust radiative transfer models con- 
strained by the observed SED and submillimcter and mil- 
limeter radial intensity profiles, L673-7-IRS is a Class 
protostar with Li nt = 0.01 — 0.045 L Q , with Li nt ~ 0.04 
Lq the most likely value. Based on the definition of a 
VeLLO as a protostar embedded within a dense core with 
L int < 0.1 L Q (Di Francesco et al. 2007), L673-7-IRS 
qualifies as a VeLLO. It is possibly the lowest luminosity 
VeLLO yet studied, although such a claim is difficult to 
prove with any certainty given the factor of ~ 2 uncer- 
tainty in Li n t for all VeLLOs studied to date (Young et 
al. 2004; Dunham et al. 2006; Bourke et al. 2006; Lee 
et al. 2009). Given that the dust in the envelope will 
reach its new equilibrium temperature very quickly fol- 
lowing a change in Li nt (<< 100 yr; Draine & Anderson 
1985; Lee 2007; Dunham et al. 2010), this is a very good 
measurement of the current L.;„ t . 

The properties of the molecular outflow, on the other 
hand, tell a different story. From the outflow, we calcu- 
late an accretion luminosity of L acc > 0.36 Lq. This is a 
very strong lower limit for all of the reasons described in 
$7] This difference between the above values of Lj n < and 
Lacc by a factor of at least 8 can be resolved by invoking 
non-steady mass accretion. The current mass accretion 
rate must be much lower than the time- averaged rate 
over the lifetime of the outflow. 

This is not a result unique to this source. Two other 
embedded protostars have been studied in the same man- 
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ner and shown to require lower current mass accretion 
rates than the time-averaged rates inferred by their out- 
flows: IRAM04191 {L acc /L mt ~ 25; Andre et al. 1999; 
Dunham et al. 2006) and L1251A-IRS3 (L acc /L mt > 4; 
Lee et al. 2010). There is growing evidence that such ev- 
idence for non-steady mass accretion is the norm rather 
than the exception, although more protostars driving 
outflows must be studied in this manner to evaluate 
whether such a result is truly universal. Such a uni- 
versal result would support the notion detailed in Sjl] 
that accretion is episodic and that most time is spent 
at low accretion rates. However, we caution that our 
results for L673-7 do not prove episodic accretion, they 
merely indicate non-steady mass accretion. Models of 
collapsing cores that are not singular but instead have 
flattened density profiles at small radii reminiscent of 
Bonnor-Ebert spheres (Bonnor 1956; Ebert 1955) fea- 
ture mass accretion rates that rise at early times and 
then decline at late times (e.g., Foster & Chevalier 1993; 
Henriksen, Andre, & Bontemps 1997; Vorobyov & Basu 
2005b). Such models could also explain a lower current 
accretion rate compared to the time-averaged rate, al- 
though a decline in the accretion rate by an amount suf- 
ficient to give a current accretion rate nearly a factor of 
10 lower than the average rate within the required time 
(6.1 x 10 4 yr, the dynamical time of the outflow) repre- 
sents a faster decline than generally predicted by such 
models. Future work should investigate this in more de- 
tail with higher spatial resolution, better-sampled maps 
of the L673-7 outflow that allow direct calculations of the 
accretion rate as a function of time. 

L673-7-IRS is classified as a Class protostar based 
on both the calculated Tboi and Lboi/L smm ( fl4.2[) . The 
physical Stage corresponding to a Class SED is a 
protostar with greater than 50% of its total system mass 
still in the envelope (Andre et al. 1993), although the 
effects of geometry and accretion bursts, if present, can 
often remove the connection between observed Class and 
physical Stage (e.g., Robitaille et al. 2006; Dunham et 
al. 2010). In this case, however, the connection appears 
to hold: the outflow properties infer an accreted proto- 
stellar mass of M acc > 0.07 M , whereas M env = 0.2-4.5 
Mq. Thus, M env /M acc < 2.9 - 64.3; most of the mass of 
the L673-7 protostellar system appears to still be in the 
envelope. 

While L673-7-IRS is consistent with currently having a 
substcllar mass, it is unlikely to remain as such. The en- 
velope provides a substantial reservoir of mass that can 
still accrete onto the protostar. Assuming a typical star 
formation efficiency in the dense gas of 30% (e.g., Alves et 
al. 2007; Enoch et al. 2008; Evans et al. 2009), the enve- 
lope mass range of 0.2 — 4.5 M Q indicates that 0.06 — 1.35 
M© can still be expected to accrete onto the protostar, 
as long as the current low mass accretion rate indicated 
by the low Li nt is not an indication that the accretion 
process has somehow terminated prematurely. Asymme- 
tries in the pattern of accretion onto Class protostars 
can give them drift velocities relative to their cores that 
may cause them to drift out of the regions of highest 
density and prematurely terminate accretion despite the 
presence of a significant mass reservoir (Stamatellos et 
al. 2005). Indeed, Huard et al. (2006) noted an offset 
between the position of the protostar L1014-IRS and the 
density peak of the core in which it is embedded, suggest- 



ing a drift velocity between protostar and core compara- 
ble to that predicted by Stamatellos et al. (2005). While 
both the SHARC-II 350 /im and MAMBO 1.2 mm con- 
tinuum emission peak slightly southwest of L673-7-IRS 
(by 2.5 and 5", respectively), both offsets are compara- 
ble to the pointing uncertainties of the two continuum 
maps. There is no clear evidence for an offset between 
L673-7-IRS and its core in the current data, although 
this question should be revisited with higher-resolution 
millimeter continuum and extinction maps. 

As discussed in $2j the distance to L673-7 is quite un- 
certain. Herbig & Jones (1983) concluded that the most 
likely distance is 300 pc, and we have adopted their re- 
sult in this paper. However, as their discussion indicates, 
the distance could be up to 2 — 3 times larger. The value 
of Li„t, determined either by the 70 fxia flux or by the 
radiative transfer models, will essentially just scale as 
Li n t oc D 2 . A larger distance by factors of 2 — 3 would 
thus lead to a larger L in t by factors of 4 — 9. L673-7- 
IRS would no longer be classified as a VeLLO, but, with 
Lint ~ 0.16 — 0.36 Li©, it would still be a low-luminosity 
protostar. More importantly, the discrepancy between 
L^t and L acc would remain, and in fact would increase. 
Since M acc = M &cc T d , M acc cx Ffi ow oc D, and r d cx D, 
EquationUDgives L acc oc M &cc M acc oc Mf cc T d oc D 3 . For 
distances 2 — 3 times greater than the assumed distance 
of 300 pc, L acc increases by factors of 8 — 27, increas- 
ing the discrepancy between L acc and Li n t by factors 
of 2 — 3. While the actual values of Li n t and L acc are 
sensitive to the assumed distance, the conclusion of this 
paper that the current mass accretion rate must be much 
lower than the time- averaged rate over the lifetime of the 
outflow remains unchanged even if the distance is much 
greater than the assumed value of 300 pc. 

9. CONCLUSIONS 

We have presented new infrared, submillimeter, and 
millimeter observations of the dense core L673-7 and have 
discovered a low-luminosity, embedded protostar driving 
a molecular outflow. We summarize our main conclusions 
as: 

1. L673-7 is seen in absorption against the mid- 
infrared background in 5.8, 8, and 24 fim Spitzer 
images, allowing for a derivation of the column den- 
sity profile and enclosed mass of the core indepen- 
dent of the millimeter continuum data and thus 
independent of dust temperature. Estimates of the 
core mass from the 8 and 24 /im absorption profiles 
range from 0.2 - 4.5 M . The MAMBO 1.2 mm 
emission gives a mass of ~ 2 M Q , both from a di- 
rect calculation assuming isothermal dust and from 
the dust radiative transfer models constrained by 
the MAMBO data point. These mass estimates are 
consistent with each other given the uncertainties 
in the dust opacities at 8, 24, and 1200 /im and the 
fact that only upper and lower limits to the true 
mid-infrared absorption profiles could be derived. 

2. L673-7-IRS is a Class protostar with L- m t 
0.0 1 - 0.045 L Q , with L int ~ 0.04 L Q the most 
likely value. It is thus classified as a VeLLO, and is 
among the lowest luminosity VeLLOs yet studied. 
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3. The properties of the molecular outflow suggest a 
protostellar mass accretion rate, averaged over the 
lifetime of the outflow, of (M acc ) > 1.2 x KT 6 ^^ 
M yr" 1 , close to the canonical value of~2xl0~ 6 
M yr -1 for the inside-out collapse of a singular 
isothermal sphere (Shu, Adams, & Lizano 1987). 
Accretion at this rate over the dynamical time of 
the outflow indicates that the protostar has grown 
to at least 0.07 M© in mass. 

4. The resulting accretion luminosity from accretion 
at this rate onto a protostar with this mass is 
L aC c > 0.36 Lq. The discrepancy between Li nt 
and L acc suggests that the current accretion rate 
is much lower than the time-averaged rate over the 
lifetime of the outflow, indicating non-steady mass 
accretion. To date, two other Class protostars 
have been shown to have similar mismatches be- 
tween L acc and L int . 

5. Although L673-7-IRS is consistent with currently 
being substellar, there is a substantial mass reser- 
voir (0.2 — 4.5 M ) remaining in the envelope for 
further accretion, making it highly unlikely that 
L673-7-IRS will remain substellar. The high enve- 
lope mass relative to the accreted protostellar mass 
(M env /M acc < 2.9 - 64.3) indicates that L673-7- 
IRS is not only a Class protostar based on its 
SED, but also a Stage protostar with greater than 
50% of its total system mass still in the envelope. 

While the results for L int should be robust to addi- 
tional data (assuming the distance is correct), and the 
discrepancy between L acc and Li nt should be robust to 
both additional data and larger possible distances, fu- 
ture observational studies of L673-7 are needed on sev- 
eral fronts: (1) A larger area should be mapped in CO 
in order to discern the true extent of the outflow. (2) 
Higher-resolution interferometer CO maps are required 
to study the outflow on smaller spatial scales close to 
the protostar, probing the more recent mass accretion 
history. (3) Submillimeter/millimctcr continuum inter- 
ferometer data are needed in order to probe basic disk 
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